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In this paper we present a very general theoretical framework for addressing fermionic superfluids over the 
entire range of BCS to Bose Einstein condensation (BEC) crossover in the presence of population imbalance 
or spin polarization. Our emphasis is on providing a theory which reduces to the standard zero temperature 
mean field theories in the literature, but necessarily includes pairing fluctuation effects at non-zero temperature 
within a consistent framework. Physically, these effects are associated with the presence of pre-formed pairs (or 
a fermionic pseudogap) in the normal phase, and pair excitations of the condensate, in the superfluid phase. We 
show how this finite T theory of fermionic pair condensates bears many similarities to the condensation of point 
bosons. In the process we examine three different types of condensate: the usual breached pair or Sarma phase 
and both the one and two plane wave Larkin- Ovchinnikov, Fulde-Ferrell (LOFF) states. The last of these has 
been discussed in the literature albeit only within a Landau-Ginzburg formalism, generally valid near T c . Here 
we show how to arrive at the two plane wave LOFF state in the ground state as well as at general temperature 
T. 

PACS numbers: 03.75.Hh, 03.75.Ss, 74.20.-z 



I. INTRODUCTION 



The subject of superfluidity in ultracold trapped fermionic 
gases is an exciting field not only for 

its implications in atomic physics, but also because there are 
important implications for condensed matter systems, includ- 
ing perhaps high temperature superconductors |9, 10]. There 
are two important aspects which are particularly notable about 
these systems, from the perspective of the present paper. They 
can be tuned in various ways which are not available to nature- 
made superconductors. Thus, one can study the entire regime 
from BCS to Bose-Einstein condensation (BEC), simply by 
the application of a magnetic field in concert with a Feshbach 
resonance. Moreover, one can vary the concentrations of the 
two spin species arbitrarily fTTl in a fashion, this 

simulates the application of a magnetic Zeeman field. This 
latter tuneability has important implications for other sub- 
disciplines in physics such as dense QCD and (isospin asym- 
metric) nuclear matter fl4l fl5l flrill . Equally important is the 
fact that there is a rich collection o f exp erimental data from 
two different atomic physics groups I 111 fl2l Il3ll on 6 Li gases 
near unitarity, with which one can compare various theoretical 
results. 

The goal of this paper is to present an overarching theo- 
retical framework for dealing with fermionic gases over the 
entire range of BCS to BEC crossover and the entire range of 
temperature, as well as the entire range of population imbal- 
ance. The foundations of this theory lie with the initial ob- 
servation of Eagles 11711 and of Leggett 1 18] that the BCS-like 
wavefunction has a much greater generality than was origi- 
nally recognized at the time of its proposal. It is capabable 
of describing both BCS and BEC like systems, providing the 
pairing attraction is tuneable from arbitrarily weak (BCS) to 
arbitrarily strong (BEC) and one self consistently solves for 
the fermionic chemical potential. This mean field-like ground 
state wave function is also readily generalized to include pop- 



ulation balance. Indeed, there are at least three well studied 
phases I19ll20l 12111 which have been proposed to accomodate 
a difference in the population of the two spin species. What 
we want to stress is that these same mean field theories have a 
natural extension to finite temperature 1 9 ] . This extension will 
be a focus of the present paper. 

At str ictly zero temperature, there is a rather extensive lit- 
erature |22[ Q ||1 |3 IH HI on these population- 
imbalanced superfluids and superconductors, including re- 
cently the effects of the crossover from BCS to BEC. There 
have also been some studies at finite T, which are at the mean 
field level and do not include the effect of the non-condensed 
pairs we consider here. The three most studied 

phases are the "Sarma"-like 1 19] or breached pair states in 
which (as in BCS theory) the condensed pairs have zero net 
momentum, but polarization can, nevertheless, be introduced. 
At T = this state appears to be stable in the deep BEC 
regime. Additionally, two different phases 12 ill have been pro- 
posed by Larkin and Ovchinnikov and by Fulde and Ferrell 
(LOFF) in which the condensate has a net momentum of a 
pre-determined q or of ±q. Even more elaborate crystalline 
lattices of various q; have also been contemplated 113 211 . It is 
believed that these states are more appropriate closer to the 
BCS side of resonance, although rather little has been deter- 
mined about the "two-plane wave" LOFF state near T = or 
in the presence of crossover effects. Added to this complex- 
ity is the possiblity of heterogeneous or phase separated states 
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In this paper we will present the theoretical formalism for 
the Sarma and one and two-plane wave LOFF states at zero 
and finite T as one varies from BCS to BEC. We note that be- 
cause theories of population imbalanced superfluids are (cur- 
rently, without exception) based on BCS-Leggett-type lTl8Tl 
ground states, it is important to determine their finite temper- 
ature implications within this broad class of ground states, as 
we do here. Our premise is that the effects of finite T, which 
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necessarily must be accomodated in any comparison with ex- 
periment, must be compatible with the T = formalism. In- 
deed, one of the most important effects of temperature is to 
stabilize the Sarma-like phase. In this way one finds an inter- 
mediate temperature superfluid I33 U341 135U36II . one that ex- 
ists at T 7^ 0, but not at T = 0. Our studies of the two-plane 
wave LOFF state present new results by extending the current 
literature away from the Landau-Ginzburg regime (near T c ). 

With this class of generalized mean field theories we will 
show that the effects of temperature enter in a fashion which 
is strongly reminiscent of Bose gas condensation. Except in 
the BCS limit, pairs form at higher temperatures, than the 
(transition) temperature T c at which they condense. Thus 
we have to distinguish the excitation gap [called A(T)] from 
the order parameter [called A SC (T)]. Alternatively, this 
means that there is an excitation gap (or pseudogap A 2 g ) for 
fermionic excitations even in the normal phase. It also im- 
plies that below T c there will be additional modes of excit- 
ing the condensate- via pair excitations. These latter are non- 
condensed or incoherent pairs with finite center of mass mo- 
mentum. We have found |34] that in a trapped geometry they 
are particularly important for providing a mechanism of get- 
ting polarization into the gas. 

These effects of finite T can be compared with an alter- 
native class of theories in the literature based on work by 
Nozieres and Schmitt-Rink (NSR) 13711 . This approach is 
known to lack self-consistency 1 38 ] . We stress that the finite 
temperature NSR approach was not designed to be consistent 
with the standard ground state equations. This observation has 
also been made in Ref. 13911 . Indeed, these latter authors have 
presented in considerable detail 1 4-0] a more complete finite T 
approach based on the T c calculations of Ref. 13711 . A major 
concern about this class of theories remains to be addressed. 
Because the gap and the number equations are not treated on 
an equivalent basis, it is possible that the superfluid density 
will not consistently vanish at T c . At this temperature there 
has to be a precise, but delicate cancellation of paramagnetic 
and diamagnetic current contributions, as found in the present 
theory 14ll I42II . In a related fashion, pseudogap effects (as- 
sociated with the presence of non-condensed pairs) appear, 
within an NSR-based approach, in the number equation but 
not in the gap equation. 

We begin at the more physical level by stressing the 
analogy between condensation in this composite boson or 
fermionic superfluid and condensation in a gas of ideal point 
bosons. Our theory treats self-consistently two-particle and 
one-particle Green's functions on an equal footing. Because 
the physics is so simple and clear, we can fairly readily an- 
ticipate the form of the central equations of this BCS-BEC 
generalization of BCS theory. It is important to stress, how- 
ever, that these equations can be derived more rigorously from 
a truncated series of equations of motion for the appropriate 
Green's functions |43[. 

There are three principle equations which govern Bose con- 
densation: the vanishing of the bosonic chemical potential at 
all T < T c is the first. Throughout this paper we will refer 
to this condition as the "BEC condition". It is related to the 
usual Thouless criterion, but the latter is generally associated 



only with the temperature T c . The second equation is the bo- 
son number equation. All "bosons" must be accounted for 
as either condensed or non-condensed. The third equation is 
the number of non-condensed "bosons", which are created by 
thermal excitations. This is determined simply by inserting 
the known excitation spectrum of the excited pairs or bosons, 
into the Bose distribution function. With this equation, and the 
first equation, one can then deduce the number of condensed 
bosons. 

These three central equations for bosons are indicated in Ta- 
ble I, on the far right, for true point bosons, and in the second 
column for the composite bosons which appear in fermionic 
superfluids. For these composite bosons the quantity which 
provides a measure of the "number" of bosons (N) is given 
by A 2 (T) (up to a constant coefficient, Z). This is reason- 
ably easy to see. In the fermionic regime, when the fermionic 
chemical potential is positive, A 2 (T) represents the square of 
the excitation gap. This is the energy which must be supplied 
to break apart the pairs. Thus, A 2 (T), in some sense then, 
reflects the number of pairs. How does one quantitatively 
establish the appropriate "boson number" for the fermionic 
case? This is determined via the self consistent gap equa- 
tion for A(T), which, in turn, is determined using the first 
condition: that the pair chemical potential is zero at and be- 
low T c . How does one compute the number of excited pairs? 
Once the gap equation is interpreted in terms of the appropri- 
ate non-condensed pair propagator, then one knows the related 
excitation spectrum fl q of this propagator. 

The quantity Z which appears in the last equation of the Ta- 
ble (for the composite bosons) gives the relation between the 
gap associated with non condensed pairs (A 2 S ) and the num- 
ber of pairs (J] b(Q, q )), It can be readily calculated in this 
theory; once one has the non-condensed pair propagator, Z 
appears as the inverse residue. (Deep in the BEC regime, Z is 
relatively simple to compute, for here the boson number den- 
sity approaches the asymptote n/2, where n is the fermion 
density). More precisely, the total number of bosons in the 
present case has to be determined self-consistently through 
the gap equation. It also involves the fermion number equa- 
tion through the related fermionic chemical potential. In this 
last context, it should be stressed that there is one important 
aspect of the fermionic superfluids, which is not apparent in 
Table I. For BCS-BEC crossover, it is essential to derive the 
self consistent equation for the fermionic chemical potential; 
in this problem the fermions are the fundamental statistical 
entity. This can be readily accomplished within the same 
framework used to arrive at the gap equation. The vanish- 
ing of the pair chemical potential is associated with a partic- 
ular choice for the pair propagator involving dressed Green's 
functions. These, in turn, determine the fermionic chemical 
potential through the fermion number equation. In the next 
two sections we turn to the gap and number equations, and 
show through a Green's function formulation, how strongly 
these two equations are inter-connected. 

The rest of the paper is organized as follows. We conclude 
this section with a summary of the central equations associ- 
ated with our T-matrix scheme. In Sec. [H] we present a 
mean-field theory for the Sarma and one-plane-wave LOFF 
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Composite bosons 


Point bosons 


Pair chemical potential 


f^pair — 0, T ^ T c 

Leads to BCS gap equation for A(T) 


Ms = 0, T < T c 


Total "number" of pairs 


A 2 (T) = A 2 sc (T)+A 2 pg (T) 


N = N + N T 


Noncondensed pairs 


zA 2 pg = E q ^o Kn*) 





Table I: BCS theory by way of BEC analogy. Here we compare condensation in composite and point bosons; \ib is the bosonic chemical 
potential, No is the number of condensed and Nt is the number of non-condensed bosons. We define \i va ir as the chemical potential for the 
non-condensed pairs. Here A(T) is the total fermionic gap which contains contributions from the non-condensed (Ap g ) and condensed terms 
(A 2 C ). In the strict BCS limit A pg — 0, so that the order parameter and gap are identical. 



state for T < T c . Due to the complexity of the two-plane- 
wave LOFF state, we dedicate an entire section fSeclIIIb to its 
mean-field treatment. In Sec. II VI we present a generalization 
of our T-matrix formalism to include pairing fluctuation ef- 
fects. Section|V]recapitulates our simple physical picture and 
Section fvTl presents our conclusions. Additional, more tech- 
nical details are given in two Appendices for the one and two 
plane wave LOFF states. 



A. Central Equations of T-matrix scheme 

To make contact with the general class of mean field the- 
ories (including the Sarma and LOFF states), we introduce 
a T-matrix approximation. This means that we consider 
the coupled equations between the particles (with propagator 
G) and the pairs [with propagator t(P)] and drop all higher 
terms. This theory does not include direct "boson-boson" 
interactions, although the pairs do interact indirectly via the 
fermions, in an averaged or mean field sense. Throughout 
Sections ITI1 and Hill of this paper we will be showing that the 
BEC condition noted in the previous section will give the 
same gap equation we find using standard techniques, such 
as Bogoliubov diagonalization applied to the linearized mean 
field Hamiltonian. Here, for all T < T c , the BEC condition is 
interpreted as requiring that the pair chemical potential /i pa ; r 
associated with the non-condensed pairs vanish. In SectionslTTI 
and[H1]we will address only the first line of Table I. The sec- 
ond two lines, or sets of equations will be discussed in Section 

EvJ 

Within a T-matrix scheme, the pair propagator is given by 

r 1 (P) = u- 1 + x(P) (i) 

where \ 1S the spin symmetrized pair susceptibility, and U < 
is the pairing interaction strength. The function x(T) is, in 
many ways, the most fundamental quantity we introduce in 
this paper. It provides the basis for obtaining well known (as 
well as new) results of the zero temperature theory. More- 
over, it provides the basis for arriving at a finite temperature 
description, which appears in Section lTVl The introduction of 
spin symmetrization is only important for the case of popu- 
lation imbalance. In earlier literature |9], this complexity did 



not arise. We will show that one obtains consistent answers 
between T-matrix based approaches and standard mean field 
theories, provided the components of the pair susceptibility in 
the presence of population imbalance, defined by 

x(P) = l[xn( p ) + xii(P)] (2) 

is given by the product of one dressed and one bare Green's 
function 

Xu(P) = J^GoriP - K)Gi(K) (3a) 

K 

Xtf(P) = ^G oi (P-^)G T (^) (3b) 

K 

where P = (ifi;,p), and G and Go are the full and bare 
Green's functions respectively. We will discuss G in more 
detail on a case by case basis. Here Gq~* (if) = iuj n — £k,o-, 
£k,<r = ek — [i a , ek = h 2 k 2 /2m is the kinetic energy of 
fermions, and \i a is the fermionic chemical potential for spin 
a =|, j. Throughout this paper, we take H = 1, kg = 1. and 
use the four-vector notation K = (io;„,k), P = (iS!;,q), 

J2 K = T S„Sk> etc ' where = (2n + l)7rT and 
fi; = 11-kT are the standard odd and even Matsubara frequen- 
cies ll44ll (where n and / are intergers). 

For the mean field discussions in Sections ||l]and|ni|of this 
paper we will not be considering general values of P but only 
zero frequency limits with special values of p associated with 
Sarma (p = 0) and LOFF (p = q ^ 0) states. However, 
when we include the contribution of non-condensed pairs (or 
pseudogap effects) the general values of P become important. 
For the Sarma phase we have the BEC condition 

tr 1 (Q) = = U- 1 + X (0) (4) 

More generally, for LOFF-like states we have the BEC condi- 
tion at finite q: 

^ 1 +X(0,q)=0 (5) 

We will discuss, in considerable detail, the nature of the mean 
field self energy which appears in the full Green's function 
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G a {K). Not only does this determine the gap equation but 
it also leads to the number equations which can be written in 
terms of Green's functions as n a = J2k G a {K). 

The number and gap equations then provide the underlying 
basis for the mean field approach. And we will see that the 
same propagator for non-condensed pairs [t(P)] enters into 
the beyond-mean-field corrections. It is important to stress 
that when we refer to "mean field" based approaches in sec- 
tions [n| and |In] of this paper, we will not be distinguishing 
between the order parameter A sc and the excitation gap A. 
Subsequently we will show that this distinction is actually an 
important one in all but the BCS limit. Specifically, we note 
that the expressions we present in Sections ITT1 and Hill within 
our Green's function-based formulation, are more generally 
valid below T c , but for the excitation gap A, not for the order 
parameter. 



function. The coupling constant U can be replaced in favor 
of the dimensionless parameter, 1/kpa, via the relationship 
m/{A-Ka) = l/U + X)k(^ £fe ) _1 ' wnere a is the two-body 
s-wave scattering length, and kp is the noninteracting Fermi 
wave vector for the same total number density in the absence 
of population imbalance. Therefore the gap equation can be 
rewriten as 



771 v - 

47TO _ ^ 



1 - 2f(E k ) 1 
2E k 2e k 



(9) 



This is a familiar gap equation which has appeared many times 
in the literature. The mean field number equations can be 
readily deduced 

n CT = ^2if( E ^)ul + f(E krf )vl], (10) 



II. GAP AND NUMBER EQUATIONS OF SARMA AND 
ONE PLANE WAVE LOFF PHASES 

A. Mean Field Sarma State 

We now want to study the Sarma or breached pair state. 
On the fermionic side of resonance, the dispersion relation of 
one of the two quasi-particle bands {E k ^ and Ek,i, defined 
below) has two zero crossings, at momenta ki (> 0) and k 2 
(> fci). This is associated with gapless excitations. In a mo- 
mentum space representation, pairing is confined to k < k\ 
and k > k 2 - This pairing is "breached" by a normal compo- 
nent in the intermediate region k% < k < k%. It is this normal 
component, then, which carries the bulk of the polarization. 

If we consider pairing between k and — k states, the mean 
field Hamitonian can be rewritten as 

h mf = J2{&A,^A+&,l c -u,iC-Kl 

k 

+ AcL kti 4,T + Ac M c -k,l } . (6) 

where we have chosen A* = A to be real. Using standard 
Bogoliubov diagonalization techniques, we readily arrive at 
the mean field gap equation via the self-consistency condition 



A = E U i C k,t C -k,l) 



(7) 



A" 



which can be written as 



i 

u 



l-f(E kl )-f(E^) 



2E k 



1 

U 



E 1 



V(E k ) 



2E k 



(8) 



Here /i = (/if + Hi)/ 2 and h = (p 

S,k — £fe 

[f(x + h) + f(x — h)]/2, where f(x) is the Fermi distribution 



= —ii -f rj k anu £j k i = h + E k , where 
/i. In addition, we define the average f(x) = 



where a = —a and the coherence factors v\ — (1 ± 
E] si )/2. Equivalently, they can be rewritten as 



E k ' 



(11a) 



Y^[f(E k -h)-f{E k + h)], (lib) 



where 77 = n^+rii is the total atomic density, Sn = n^ — n^ > 
is the number difference and p — Sn/n is the polarization. 



B. T-matrix Approach Below T c : Sarma State 

We next show that the same results can be obtained from a 
T-matrix based approach. This discussion provides the link 
between the first line of Table I and the previous subsection. 
The one particle Green's function for particles with spin a is 

G-\K) - G^{K)-^{K) 

= iu n - £ feCT - E<r(K) , (12) 

The self-energy Y, a can be shown to be of the BCS-like form 

A 2 



(13) 



We will see later in Section llVl how this form for the self en- 
ergy very naturally arises (below T c ) in a T-matrix approach. 
Thus 



G-\K) =kj-£ kcr - 



A 2 



7W + S,kS 



(14) 



Then, using the coherence factors we defined earlier, the 
Green's functions become 



G a (K) 



ilo - E krT iu> + E ks 



(15) 
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Now we are in position to calculate the pair susceptibility 
at P — for the Sarma phase based on Eqs. (0 and 



x(o) = xti(0) = xit(0) 
= ~^77 



i 



K 



(iuj n - E kl )(iuj n + E^) 



(16) 



Substituting this expression into our BEC condition Eq. ©, 
we obtain the same gap equation (|8), after carrying out the 
Matsubara summation. 

In terms of Green's functions, we readily arrive at the num- 
ber equations: n a = J2k G a {K), which reduce to the num- 
ber equations (II Oi we found earlier. 



C. Mean Field Theory of One Plane Wave LOFF State 

If we now consider condensates in which momentum k 
pairs with — k + q, for, as yet undermined q, the mean field 
Hamiltonian can be rewritten as 



H 



MF 



{&>T C k,T Ck >T + £k-q,i C -k+q.l C -k+q,l 

k 

+ Ac ik +q ,i c k, T + Ack, T c_k +q ,i} d7) 



Upon Bogoliubov diagonalization, the self-consistency condi- 
tion 



(18) 



K 



readily leads to the mean-field gap equation 
U ^ 



2Ekq 



(19) 



with Ei. 



E kq ~h+ (e fc - e k - q )/% E ltl = E kq + h 



(e k - e fc _ g )/2, E kq = ^ + A 2 , £ kq = (& + &_,)/2. 
The regularized gap equation is 



m \ - 

A-KCL 4"^ 



1- f{E^)- f{E n ) 
2E kq 



(20) 



Finally, the mean field number equations for n a , are given by 
n a = Y,{f( E ^>l q + f(-Ei.a)v 2 kq }- (2D 



Here the coherence factors are u% q = | (l + JrM and v kq 



2\ E k g t 

There must be another equation which governs q. This can 
be obtained by minimizing the thermodynamical potential or 
free energy with respect to q. Equivalently, we will derive q 
from the T-matrix method described below. 



D. T-matrix Approach Below T c : One plane wave LOFF state 

We now use the same T-matrix based approach (as we did 
for the Sarma case, and as indicated by Table I) to make con- 
tact with the one plane wave LOFF state. We take the self 
energy of the form 



= -A 2 G a9 (-K) 



A 2 



IU) + Ck-q.ff ' 



so that 



G-\K) = kj-i 



A 2 



k,<r 



k— q,cr 



Then we have 



G T (K) 
Gi(K) 



iu) n - Ei^ iuj n + Ei }i 



iuj n - E 2 ,i 



iuj n + E- 



(22) 



Here T 2 ,t = E kq ~h~(e k ~ e k - q )/2 and E 2 .\ = E kq + h + 
{(-k — efe-q)/2. Note if k — > — k + q, then E 2 .\ — > Pij and 
E 2t i — > The number equations are n a = J^k G a (K), 

which yields the same answer ("Eg. 12 lb we found in the mean 
field approach. For T < T c , using the BEC condition 
U^ 1 + x(0,q) = we thus arrive at the gap equation we 
found earlier in Eq. dl9l . 

We next investigate the pair susceptibility and its extremal 
value at Q — (0, q). We define the contributions to the pah- 
susceptibility in the presence of population imbalance as in 
Eqs.(Q]i-(|3- For the one plane wave LOFF state, the BEC con- 
dition is that t(Q, p) diverges at a non-zero momentum p = q. 
Thus, at this momentum, the quantity x(0, p) should reach a 
maximum at p = q. We determine q by requiring 



Bx(0,p) 



dp 



o. 



(23) 



Where x(P) is explicitly shown in Eq. (IA1> . This yields 

0X(O,P)| 



= 



dp l P = 
A 2 ^ | 2 



l-p-)-[f(Ei,t) + f(Ei,i)}p- 

&kq ' &kq 



(k-|) lf(Ei A -)-f(Ei d )} 



(24) 



It is important to stress, as we show in Appendix [A] that 
this extremal condition on x(P) is equivalent to the condition 
that the net current is identically zero in this situation. 

j = Y,k[G r (K) + G l (K)} = (25) 

K 

This key observation shows that the present way of computing 
q directly from the pair susceptibility is consistent with the 
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counterparts in the literature, which are based on the vanishing 
of the net current in equilibrium. 

At the mean field level we, then, have four unknowns A, q, 
ji and h, and four equations: two number equations, the gap 
equation, and the condition on the vanishing of the first order 
derivative. 



III. MEAN FIELD THEORY OF TWO PLANE WAVE LOFF 

STATE 

The two plane wave LOFF state has not been studied in 
much detail away from T c and the associated tri-critical point 
ll32ll . Here we address this state at more general temperatures, 
starting with the very natural mean field Hamiltonian in the 
same spirit as in the previous two cases. The Hamiltonian is 
given by 



H 



MF 



E 

k 



^ A ( C kT C -k+ql + CkiC_ k+qT ) 



v/2 



A ( C kT C ^k-qi +C kl C -k-qT 



)} (26) 



To solve for the value of the gap we generalize the Bogoli- 
ubov diagonalization procedure. If one rewrites the effective 
Hamiltonian in the basis set = (ekf, c_ k _ q ^, cL^qj^) 

and — (ck|, cl k |, c -k+q|)' tne Hamiltonian can be 
written in two 3x3 matrices. Here we show the matrix for 



B 



Ht 



64 



7I A 

~£k+qi 




7I A 




-6 



(27) 



-qi / 



We can equivalently write the | component of the 3x3 

matrix in terms of the basis H'J = (ck+ q | , Qc-qt , c !ki ) The 
mean field gap equation is derived from the averages obtained 
from a symmetric combination of the two basis sets as 



k,l) 



K 



which is equal to 



K 



K 



K 



U (c k+q . 



T c -fc,l/ 



(28a) 



(28b) 



We choose this symmetrized form to write down the gap equa- 
tion, which reflects the underlying symmetry of the Hamilto- 
nian. 

The resulting gap equation can be shown to be of the form 



1 

U 



-E 
2 ^ 

k 



f( El1 ) + f(E H 
±[f(E 21 )+f(E 2l )] 



E x + £ k . 



kq 



(Ei — E 2 )(Ei 

E 2 + £ kq 



E 3 ) 



(E 2 



lf(E 3 i) + f(E 3l )} 



-Et)(E 2 - 

E3 + £kq 



E 3 ) 



{E 3 -Ex)(E 3 -E 2 ) 



■(29) 



Here we define E±, E 2 , and E 3 as the solutions to the cubic 
equation which is related to the determinant of the 3x3 eigen- 
value equation 



-e. 



kq-i 



E kq- 
(k 2 - 



kq 



5-£ 



kq 



E: 



kq 



(—)' 

\ m J 



= 0. 
(30) 



\q 2 )/2m - fi and E kq = J £ 2 q 



where Ekq 

Then E n = Ej - (q 2 /4m) - h and E n = Ej - (q 2 /Am) + h 
forj = 1,2,3. 

The mean field number equations are 



+ a-vlq-l 2 kq)f(E 3 ,) 



where we define 



vtq 



llq 



(E 1+£kq ) 2 -(^) 2 
(E 2 — Ei ) (E 3 — Ei ) 

(E 2 +£kq) 2 ~(^) 2 

(Ei — E 2 )(E 3 — E 2 ) 



(31) 



(32) 



A. Green's Function Approach Below T c : Two plane wave 
LOFF state 



For the counterpart mean field theory, approached from a 
T-matrix scheme, we define the self energy as 



Z a (K) = -^[a 2 G 0s (-K + Q) + A 2 G 0S (-K - Q) 



A 2 



A 2 



£k— qcr 



£k+qer 



(33) 



HereQ = (0,q). 

The Green's function is then given by 



G a (K) = 



1 



GoAK)-Z a (K) 

(juj n + £k+qg)(^n + £k-qtr) 

(ioj n - Ei„)(iu)„ - E 2(T )(iuj n - E 3cr ) 



(34) 



Tfco 



1 - Vkq ~ llq 



iuj n - Ei a iu„ - E- 



2a 



E : 



Again, following Eqs. Q-Q we can write out the form of 
the pair susceptibility. For the two plane wave LOFF state, 
this equation is presented as Eq. dBll in Appendix [B] This 
quantity, in turn, enters the gap equation, given by the BEC 
condition, 1 + Ux(0, q) = 0. We may write this gap equation 
in compact form as 



—E 



£kqer 



2 ^ (iuj n - Ei a )(iuj n - E 2a )(iuj n - E 3a ) 



(35) 



This will, in turn, reduce to the gap equation ( 1291 we deduced 
directly from the 3x3 matrix analysis. 
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Similarly, in our Green's function formalism, we have 
n a = ^2 K G a (K), which reduces to the number equations 
we found in Eq. J3 li . 

Next we determine the momentum p that maximizes x(P) 
when fl — 0, i.e., we need to find a solution to dx(°>p) = 0. 

dp 

Since the Green's function is symmetric under k — > — k, it can 



be shown that p = is a solution, which corresponds to the 
Sarma phase. Here, however, we are interested in a LOFF-like 
state p = ±q. where both signs contribute in a symmetric 
fashion. Thus we choose p = q, and, thereby, arrive at the 
defining equation for the net momentum of the pairs. 



p / k - q \ f iuj n +gk+qff 1 "I 

7^ V m J\ lUJn + £k-q<r (lOJn - E\„){<bW n - E 2a )(iuj n - E 3a ) ) 

p/ k-q u 2(±Z) [/(&- qT ) + J(g k - qi ) - 2] 

?\ m )\ {El +S kq -^)(E 2 + S kq -^)(E 3 +S kq -^) 

f(E^) + f(E u ) E 1+ £ kq + ^ f(E 2] ) + f(E 2i ) E 2 +£ kq + ±* 
( El - E 2 )( El - E 3 ) Ei + £ kq - ^ (E 2 - Ei)(E 2 - Eg) E 2 + £ kq - ^ 

f(E^) + f(E 3l ) E 3 +8 kq + ^ ^ 
(E 3 -Ei)(E 3 -Ei) E 3 + 8 kq ~^y 



This determines the magnitude q; for definiteness we take the 
direction of q as along the z-axis. 

The case of the single plane wave LOFF state should be 
contrasted. There we showed that there was an intimate re- 
lation between the condition that there be no net equilibrium 
current and the extremal requirement on x(-P). For the two 
plane wave LOFF state the current can be shown to be identi- 
cally zero. Rather, the only condition one has to determine q 
is the vanishing of the first derivative of the pair susceptibil- 
ity. A similar condition was imposed in the original Larkin- 
Ovchinnikov (LO) paper 112 ill . 



(a) (b) 




= (1/2) [ > » < + > ■ < ] 



Figure 1 : Diagrams for Sarma states or alternatively one plane wave 
LOFF states including (a) Bare Green's function Go- (b) Vertices 
from pairing. The momenta are (k j, — k ]) and (k "f, — k j), while 
for single plane-wave LOFF states the momenta are (k j, — k + q t) 
and (k f, — k + q |). (c) Green's function G. (A) (ck \ c_k|) 
for Sarma states and (c k | c -k+ q |) for single plane-wave LOFF 
states. Here the thin and thick lines represent bare and full Green's 
functions, respectively. 



B. Diagrammatic Interpretation of Green's Function 
Approaches: Comparison with Larkin-Ovchinnikov 

We now want to compare the present approach for the two 
plane wave LOFF state with that proposed in the original LO 
paper l2lll . We will do this comparison within a diagram- 
matic framework. It is useful first to illustrate the diagram- 
matic scheme by referring to the simpler Sarma and one plane 
wave LOFF states. These two states have a rather similar dia- 
grammatic formulation. We write these diagrams in a consol- 
idated form in Fig. [2 The first line indicates the bare quanti- 
ties, that is, the bare Green's function Go in part (a) and the 
quantity Go^(K)AGo y i (— K) in part (b). This is written for 
the Sarma state and readily generalized to the one plane wave 
LOFF phase. 

The second line indicates the (diagonal) dressed Green's 
function, while the third line shows effectively the Gor'kov F 

Go ki -k+qt k( -k+q| M -k-qi kl -k-qf 

(d) 

* = > + > * < * > + > » ~ < * > + 



(e) 




Figure 2: Diagrams representing the original Ovchinnikov-Larkin 
theory for (a) Non-interacting Green's function Go. (b) Vertices from 
pairing with q. (c) Vertices from pairing with q. (d) Green's func- 
tion G. (e) (c k T c-k+qi)- 
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w W + + to 

Go kt -k+qt k( -k+qt kf -k-ql k( -k-qf 

(d) 

+ + _ _ 
S = > + > • — < — • > + > « < — • > + ... 

+ + — — 

= > • < ■ > + > • — < ■ > 

(e) 



= (1/2) [ > ■ < + > • < ] 

Figure 3: Diagrams associated with the 3x3 Bogoliubov diagonaliza- 
tion theory for (a) Non-interacting Green's function Go- (b) Vertices 
from pairing with q. (c) Vertices from pairing with — q. (d) Green's 
function G. (e) (c k T c_ k + q j.)- 

function. It should be noted for Fig.^1 that this latter "anoma- 
lous" or off-diagonal Green's function is explicitly seen to de- 
pend on a symmetrized sum of the product G and Go- It is this 
combination which we have seen appear in the pair suscepti- 
bility x(-P). Indeed, throughout this paper we have found that 
GGq is essential for arriving at the standard mean field theo- 



retic approach. We can, thus, conclude from the last line in 
the figure that the F function is given by a spin symmetrized 
combination of AGGo- 

In Fig.|2]we show the two plane wave LOFF diagrams orig- 
inally proposed by LO. In order to discuss their implications, 
we present the diagrams associated with the mean field matrix 
scheme of the previous section which are plotted in Fig[3] As 
in Figure^ one can see from the last line in part e of Fig[3]that 
the off-diagonal Green's function depends on a symmetrized 
sum of G and Go- One can also see that there are similari- 
ties as well as differences in these latter two approaches. In 
the approach of Section Ull Al just as in LO we restrict our 
calculations to diagrams that do not contain propagators with 
momentum beyond ±k and ±k ± q. The differences are also 
apparent. In the LO paper, the (diagonal component of the) 
Green's function diagrams were summed up to second order 
in A. In the gap equation LO dropped diagrams containing 
propagators with higher momenta and truncated the series at 
third order in A. 

If one were to follow the original LO scheme, but sum the 
entire series, then we will arrive at a slightly modified gap 
equation. We would find instead the gap is determined by 



1 

U 



]T {[G 01 (-K + Q)G l (K) + G 0L (-K + Q)G ] {K)] - \[G^(-K + Q)g{ 1) (K) + G 0L (-K + Q)G\ 1 \k)}} 

K 

£{[/(£ 1T ) + /(£u)] El + £kq 1 — - 1 - E2 + £kq 

k 

+ [/(£ 3T ) + f{E H )\ 



(Ei — E2){E\ — E3) ^ 2 ^ ^ 2 ^ (E 2 — Ei)(E 2 — E3) 



E z +£. 



kq 



(Es — Ei)(Es — E2) 



2E 



(i) 

kq 



(37) 



and thus the pair susceptibility derived from this method is 

X(P) 

= \ E { [ G °tOP - K)GdK) + G oi (P - K)G r (K) 



K 



- [g ot (P - K)G{ 1] (K) + G ol (P - K)G\ 1] (K)] }. 

(38) 

We could similarly use this modified pair susceptibility to ar- 
rive at the constraint on the value of q. We discuss these con- 
tributions in Appendix 151 For the above equation we define 
the Green's function for the single plane-wave LOFF states 

as Gf{K) = (tw n +€ k _ qi )/(iw B -£?g)(iw n + Eg)and 



G^iK) = (iw„+&_ qT )/(iw n -£^)(iw n +^). Theen- 
ergy spectrum of the single plane-wave LOFF states is e[^ = 

+h, 



^(k-q- and = E, 



a 2,-[ — ^kq 



(1) 

kq 



1 

2m 



h, where E, 



(i) 

kq 



(£kq - 



2m / 



A? 



and we define Af = 



5A 2 



In summary, to make contact with the original results of 
LO, one must subtract a second symmetrized term, which rep- 
resents pairing with only one momentum. It would thus ap- 
pear, that relative to LO, there is an overcounting in the sum- 
mation of the two series from the diagrams shown in Fig|3] 
and associated with the 3x3 matrix or mean field representa- 
tion. At this stage it is difficult to determine which of the two 
plane wave LOFF representations is the more appropriate. It 
will be essential in future to study them both numerically in 
the presence of BCS-BEC crossover effects. 
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IV. BEYOND SIMPLE MEAN FIELD THEORY: PAIR 
FLUCTUATION EFFECTS 

A. Inclusion of the Pseudogap 

For definiteness the equations which appear in this section 
apply to the Sarma-like phase. We can readily generalize to 
include the two different LOFF states. 

This diagrammatic representation of our T-matrix scheme 
is shown in Figs.|4]and|5] The first of these indicates the prop- 
agator for non-condensed pairs which we refer to as t pg , and 
the second of these the total self energy. One can see through- 
out the combination of one dressed and one bare Green's func- 
tion, as represented by the thick and thin lines. The self energy 
consists of two contributions from the non-condensed pairs or 
pseudogap (pg) and from the condensate (sc). There are, anal- 
ogously, two contributions in the full T-matrix 



tpg(P) = 



u 



1 + u x (P) : 



t sc (P) = —f5(P) : 



(39) 

T^O, (40) 
(41) 



where we write A sc = — U Ek( c - k l Ck T)- Similarly, we have 
for the fermion self energy 



H a {K) = K C (K) + TP a a {K) = £ t(P)G , s (P ~ K) 



(42) 



We can see at once that 

K C (K) = Y, t sc(P)G , s (P -K) = -G Q , S {-K)A 2 S 
p 

The vanishing of the pair chemical potential implies that 



(43) 



t- 1 (0) = U- 1 +x(0)=0, T<T C 



(44) 



Moreover, a vanishing chemical potential means that t pg (P) 
is strongly peaked around P = 0. Thus, we may approximate 
l45ll Eq. E3 to yield 



where 



A 2 (T) = A 2 C (T) + A 2 (T) . 



Importantly, we are led to identify the quantity A 

P#0 



I'U 



(45) 



(46) 



(47) 



Note that in the normal state (where [i pa ir is non-zero) Eq. 
(145 \ is no longer a good approximation. 

We now have a closed set of equations for addressing the 
ordered phase. We can similarly extend this approach to tem- 
peratures somewhat above T c , by self consistently including 
a non-zero pair chemical potential. This is a necessary step 



+ + 



+ 



+ 




Figure 4: T-matrix diagrams. Note that there is one dressed and 
full Green's function. This represents the propagator for the non- 
condensed pairs. 



E = 



Figure 5: Self energy diagram for the present T-matrix scheme, in- 
dicating the condensed (S se ) and non-condensed (E P9 ) pairs. 



in addressing a trap as well l34Tl . Additionally, the propaga- 
tor for non-condensed pairs can now be quantified, using the 
self consistently determined pair susceptibility. At small four- 
vector P, we may expand the inverse of the T-matrix, after 
analytical continuation (ifli — > fl + i0 + ), to obtain 



ai n 2 + z(n 



2M* 



f-Lpai 



(48) 



where the imaginary part Tp — > rapidly as p — > below T c . 
Because we are interested in the moderate and strong coupling 
cases, we drop the aifl 2 term in Eq. d48l >. and hence 



t P9 (P) = 



where we associate 



Z 



2M* 



(49) 



(50) 



This establishes a quadratic dispersion and defines the effec- 
tive pair mass, M*. 

Finally, one can rewrite Eq. 1471 as 



Al g (T) = Z-^ b (n p ). 



(51) 



where b(x) is the Bose distribution function. Analytical ex- 
pressions for this mass are possible via a small p expansion of 
X- In this way we find 



x(P) - x(o) 

The coefficients are 

Z = 

and 

1 



z n 



2NP 



£1=0, p=0 



1 d 2 x 



2M* 



6Z dp 2 



(52) 



(53) 



(54) 



n=o,p=o 
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B. Pair Dispersion in Sarma State and Physical Consequences 
at Finite T 

The theoretical framework in the previous section can now 
be implemented to address the physics at finite T. We do this 
here only for the Sarma state, thereby expanding on earlier 
papers 13 3l 13411 . The counterpart results for the two LOFF 
phases are presented in the two appendices. 

We begin by determining the form of the pair dispersion. 
The pair susceptibility is given by 



E 



u2 /(£fc)+/fo,-fc)-l 

k ifl - £ p -k - Ek 



(55) 



Inserting this form into the T-matrix [Eq. Q] we readily find 
that the coefficent Z is 



Z 



dx 

an 
i 

2A 2 



n=o,p=o 

-2^/(6) 



(56) 



Here n = ^ k a n a (k) is the total density. 

Next we want to calculate £l p (defined above) and its deriva- 
tives. 



1-1, 

— t fi=0 



where ((x) is the Riemann zeta function. When T > T c , 

Physically, these non-condensed pairs, which appear at fi- 
nite T, have been shown to have important physical conse- 
quences. In the homogeneous situation, we have found that 
finite T stabilizes the superfluid state leading to an "intermedi- 
ate temperature superfluid" j^lESl- In addition, the presence 
of a finite excitation gap at T c [A pg (T c ) ^ 0] makes the be- 
havior of the superfluid transition temperature more complex 
than its mean field counterpart, T C MF '. We frequently find a 
double-valued structure |33], with superfluidity existing only 
for temperature intermediate between the two T c 's. 

In the trapped situation, we have found 1 34] that these pairs 
enter explicitly as the mechanism for carrying polarization 
within the Sarma phase. At low T there is very little polar- 
ization carried by the superfluid core; instead, the region of 
the trap where A sc = 0, but A / - which can be called 
"the mixed normal region"- carries the bulk of the polariza- 
tion, much as observed experimentally H 1 311 . 



V. PHYSICAL PICTURE: FERMIONIC PAIRING IN 
ANALOGY TO BEC 

We now return to the strong analogies between this BCS- 
based or Leggett mean field theory 1 18] and Bose condensa- 
tion of point bosons, as summarized in Table I. 

We have three central equations. 

1 . The pair chemical potential must vanish at and below T c 



1 m x - 

~Z\ 47m _ ^ 



1 



1 -/(**)- /&>-*) 



,m P - k )~i(E k ) 



Since £ fc = 
3/m. Then 

1 



Ek — £, P -k 



k 2 /2m—fj,, we have Vfc£fc 



id 2 n r 



(57) 



k/m and V|£fc 



2M* 



6 dp 2 I p=o 



12toZA 2 



£{/(£ 



6 



f(E k] )[ 



66 



mA 2 

m+ri(E 2 k +e k ) 



E k 

fa + 
m 

6 



3 1 



mE k A 2 
Ej 



(58) 



Here f'(x) is the derivative of f(x). 

With this dispersion, then one can compute the number of 
non-condensed pairs, n pair = ZA pg = YlpH^p)- ^e then 
have 



A 2 

pg 



(2M*T) 



3/2 



2n 2 Z 



(59) 



f-^pa 



0, (T<T C ). 



(60) 



Importantly this condition leads to the mean field gap equa- 
tions derived in Sections UTI and fiTTl These gap equations then 
provide a specific value for A(T), according to the different 
phases being contemplated. 

2. There must be a conservation of the total number of 
(composite) "bosons" in the system. For this condition, our 
central equation is Eq. ( I46l l. Here it is understood that the 
number of "bosons" is effectively represented by the param- 
eter A 2 (T). In the fermionic limit, this parameter reflects 
the number of bosons through the energy which is needed 
to create fermions, and thereby break the bosons apart. Un- 
like the point boson case, here the "total boson number" is 
temperature dependent and has to be self-consistently deter- 
mined. As expected, in the deep BEC regime, where the 
fermionic excitations are negligible, the pair density is given 
by n pair = ZA 2 « n/2. 

3. The number of non-condensed pairs is readily computed 
in terms of the pair dispersion, just as in conventional BEC. 
For this condition our central equation is Eq. J5 II . 

Then, just as in conventional BEC, the number of con- 
densed bosons (proportional to A 2 C ) is determined by the dif- 
ference between A 2 (T) and A pg (T). This, in turn, deter- 
mines the transition temperature T c as the lowest tempera- 
ture^) in the normal state at which noncondensed pairs ex- 
haust the total weight of A 2 so that A 2 g = A 2 . Solving for 
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the "transition temperature" in the absence of pseudogap ef- 
fects I29tl30tl3 ill leads to the quantity T C MF . More precisely, 
T C MF should be thought of as the temperature at which the 
excitation gap A(T) vanishes. This provides a reasonable 
estimate, for the pairing onset temperature T*, (when a sta- 
ble superfiuid phase exists). This is distinguished from the 
transition temperature. We note that T* represents a smooth 
crossover rather than an abrupt phase transition. 

It should be stressed that the dispersion relation for the non- 
condensed pairs is quadratic. This appears in the Sarma and 
both one- and two-plane wave LOFF states. For the latter 
two states, it can be seen to be closely related to the fact that 
x(0, q) must reach a maximum at q. While one will always 
find a linear dispersion in the collective mode spectrum M-1XI . 
within the present class of BCS-BEC crossover theories, the 
restriction to a T-matrix scheme means that there is no feed- 
back from the collective modes onto the pair excitation spec- 
trum. In effect, the T-matrix approximation does not incorpo- 
rate pair-pair interactions at a level needed to arrive at this ex- 
pected linear dispersion in the pair excitation spectrum. Nev- 
ertheless, because essentially all theories which address pop- 
ulation imbalance build on the simplest BCS-Leggett mean 
field theory, there is good reason to first address this level of 
approximation when including finite temperature effects. 

VI. SUMMARY 

It should be clear from the previous sections that the zero 
and finite temperature theories of population imbalanced su- 
perfluids can be consolidated into one general theory, based 
on the quantity which we call the pair susceptibility, x(P). 
This has a very specific form in the class of mean field the- 
ories currently applied to address population imbalance. The 
quantity x(P) enters into the propagator for non-condensed 
pairs Eq. J40i which is just the T-matrix, t(P). When non- 
condensed pairs are in equilibrium with a condensate, they 
must have zero chemical potential. This, in turn, yields the 
various gap equations for the Sarma and the one- and two- 
plane-wave LOFF states, provided one take a special form 
for the pair susceptibility involving one dressed and one bare 
Green's function. 

Importantly, the same pair propagator characterizes the ef- 
fective number of non-condensed pairs as seen in Eq. J47I . 
Each of these population imbalanced superfluids has an asso- 
ciated T ^ pseudogap contribution, which is summarized 
for the Sarma case in Section[|I1]and for the two LOFF cases 
in Appendices lAl and iBl This pseudogap contribution serves 
to differentiate the order parameter from the gap parameter at 
all non-zero T, and in all cases except strict BCS theory. At a 
more physical level, in the normal state there is an excitation 
gap for fermionic excitations associated with "pre-formed" 
pairs. In the superfiuid phase, there is a new form of con- 
densate excitation, not found in BCS theory and associated 
with excited pair states. As the pairing attraction becomes 
stronger, it pays to excite pairs of atoms rather than create 



single fermion excitations which cost an energy gap. These 
pseudogap effects are an essential component of BCS-BEC 
crossover theory, and they are necessary in order to smoothly 
involve from the fermionic statistics of BCS to the bosonic 
statistics of BEC. Interestingly, they are also widely observed 
in high temperature superconductors fl Uoll . 

Thus far, we have applied the present theoretical formalism 
to the Sarma state both in the homogeneous 1 33 ] and trapped 
ll34ll configurations. From the point of view of comparing 
with experiment, our trapped calculations have again under- 
lined the importance of pseudogap effects. We find that the 
bulk of the polarization is carried in the pseudogap region of 
the trap: outside the condensate but in the region where both 
spin states (and thus pairing) are present. Because we find 
that only modest polarizations are stable, it appears necessary 
in future to include LOFF-like condensates as well, although 
their stability must be demonstrated ll36il . This observation 
is also consistent with numerical calculations ism based 
on Bogoliubov-de Gennes theory. A key contribution of the 
present work is that it lays the groundwork for addressing the 
one and two plane wave LOFF phases, at general temperature. 

In summary, this paper has presented a theoretical formal- 
ism for the Sarma and one- and two-plane wave LOFF states 
at zero and finite T as one varies from BCS to BEC, in the 
presence of an arbitrary population imbalance. Our premise is 
that the effects of finite T, which necessarily must be accomo- 
dated in any comparison with experiment, must be compatible 
with the T = formalism. The zero temperature formal- 
ism we use here reduces to the standard one in the literature 
II22II27II for the Sarma and one-plane wave LOFF states. How- 
ever, our studies of the two-plane wave LOFF state present 
new results by extending the current literature away from the 
Landau-Ginzburg regime (near T c ). 
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Appendix A: ADDITIONAL RESULTS FOR THE ONE 
PLANE WAVE LOFF STATE 

The following two appendices are dedicated to present- 
ing additional details on the one- and two-plane-wave LOFF 
states. The results are presented in appendix form both for 
clarity and to avoid some of the more technical details in 
the main text. First, for the one-plane-wave LOFF state, 
we wish to explore the relationship between the zero cur- 
rent condition applied by Fulde and Ferrell [21] and the ex- 
tremal condition of x(0,q)- For the former, we determine 
momentum q by requiring the total current to be zero. From 

J = (Vx-V./JUG'^i) = J2k [k^T ( K ) + (K)] 
we get 
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j - ^k[/(B!, t K + /(-B u )^ - ^(k-q)[/(^K + /Kr)ti] 

k k 

= E [f ( x - Jr) + f[/(^.T) + /(^u)]f^ + (k-f )[/(£i, T ) - /(^u)]" 

, kg kq 

= 

I 

We can take q as the z direction, and all the above equations x(-P) = \ [XU (P) + XlT (P)\ ^ s 
need to be integrated over the angular variable 8 . 

The pair susceptibility, even for the simpler LOFF state is 
reasonably complex so we present it here for completeness. 



1 V [ 2 /(^i.T) + f&d) + /(€p-fc, T ) + /fo-fcq) - 2 a /(^-fc,T) +/(&-*,;) -/(^,T)-/(^ia) 



- + - (e* - e k - q )/2) 



(Al) 



Next, we characterize the pseudogap contributions which Here n is the total density. Note that the linear derivative van- 
enter via Ap g = Z~ x ^ p b(Q,p). Here Z and j^p- are deter- ishes, as was established in Section HTCl Thus we turn next to 



mined as follows. 



the quadratic term. 



at- 



ari 



f2=0,p=q 



2A 



I The coefficient of this second order term corresponds to the 

^2 \ n ~ y^[/fe-gT) + f(£k-qi)]j- (A2) inverse pair mass, 1/2M*. This is given by 



2M* 



i g 2 x(Q,P) 

6Z dp 2 
1 



12mZA 2 



E{ 2 /(^- q ) 



4£ fcg (k-q) 2 



[/(£?i, T ) + f(E hl )] 



36 9 ^ 2[(k-q) 2 +g 2 ](£; fc 2 9 +^) 1 



2[/(^i, T ) - f(Ei,l)] i 



:k-q) 2 -g 2 ]a g + -f(ek- q )(k-q) 



3 fi -i*L\ 

Ekq ' 



(A3) 



Appendix B: ADDITIONAL RESULTS FOR THE TWO 
PLANE WAVE LOFF STATE 

We begin with the results that relate to the simplest mean 
field-based scheme discussed in Section ITlI Al We then pro- 



ceed to the slightly different scheme based on a resummation 
of the LO diagrams. 



The pair susceptibility is given by 
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x(P) = W^ p - K ) G ^ K ) + G ^ p - K ) G ^ 



K 



(iuj n +Ck+ qCT )(iw„ + Ck- 



2 ^ ifi - iw„ - C P -k CT («w„ - £'i 5 )(iw„ - E 2 ^)(iu n - E 3S ) 



(Bl) 



From the pair susceptibility one may obtain the dispersion relation for the non-condensed pairs. This is important for intro- 
ducing pseudogap effects and treating T^O. 

The explicit forms of the coefficients Z and 1/2M* are 



in 



r k.q[/(a- qT ) + /(6c- q i)-2] 



+ 



f{E 2] ) + f{E 2i ) 



^){E 2 +£ kq 



^){E 3 +£ kq -^) 
f(E 3] ) + f(E 3l ) 



+ 



f(E^) + f(E n ) E 1+ £ kq + ^ 



+ 



{E 1 -E 2 ){E 1 -E 3 ) E 1+ £ kq 

E 3 +£kq + ^ 



kq 



(E 2 - E 1 ){E 2 - E 3 ) E 2 +£ kq -^ (E 3 - E X ){E 3 - E x ) E 3 + £ k , 



+ m I 
k^q / 



(B2) 



and 



2M* 



k 



k-qT 



+ /(a-qi) - 2 



+ 



£ kq ^)(E 2 + £ kq ±*)(E 3 + £ kq - 
1 1 



24mZ ^ I (Et 

^^Ei+£ kq -^ E 2 +£ kq 
£(k • q)(a-q + M)[/'(Ck-q T ) + /'(a-qi)] 



+ 



kq 



+ 



12 4 

— (k-q) + -(&_„ + //) 

m m 

1 



^3 + £fcn 



kq 



+ 



+ 



(E 1+ £ kq -^)(E 2 

f(E^) + f(E u ) 
(Ei — E 2 )(Ei — E 3 ) 

f(E 21 ) + f(E 2l ) 
(E 2 — Ei)(E 2 — E 3 ) 

f(E 21 ) + f(E 2l ) 
(E 2 — Ei)(E 2 — E 3 ) 



■£k q - -^){Es + £ kq - -£) 



(E 1+ £ kq -^r 

ri(E 2 +£ kq + ^) 



r-(£ 



k — q 



m / g ! < / ■ 



- 6 



E 2 



- £ kq 



k^q 

m 

k^ 

m 



A(£ k _ q + M )(2k+£ fc , + *a) 



- 6 



+ £ fcq 
#3 + + ^ 



-^3 + £fco 



kq 



(B3) 



For the LO-based pair susceptibility we have 

X LO (P) = J2^( P K )Gl(K) + G 0l (P K)G,{K)] \ ^[G 0T (P - K)Gf(K) + G oi (P - K)Gf{K)\ 

K K 
T y-V- f 1 (^n+Ck+q CT )(^«+gk-q CT ) \ 

n ' _iw ™ _ £p-k<7 - E ls ){iuj n - E 29 )(iu n - E 3rf ) ! 

2 „ k ^ - - (^n - 45 ) + 4!{ ) 
1 iUn + frc-qj 



+ 



} (B4) 
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The momentum q is determined by minimizing x(0, p), i.e., ax ^ p ) — at p = q. This gives the momentum equation 



= T yy (±zl)i 



n k,<r 



k q\ iuj n +^k+qo- 
rn 



iVn + 6k-q<r (lUn - Eig){lUJ n - E 2s ){lUJ n - E 3S ) 

1 1 1 



k 



m 



3 ) 



n + 6k-qT (i W „ - ) (lLU n + ) «"n + £k-ql ( M „ - ) (tW„ + ) 



m 



){ 



£(k-q)[/(&- qT ) + /(g k - q J-2] 

(Si + 5 fc , - ^)(^ 2 + £ feg - ^)(£? 3 + £ kq - ^) 

kq 



/(Sit) + /(Six) s x + + ^ /(£ 2T ) + /(£ 2i ) £ 2 + £ fe9 + 



(Si - S 2 )(Si - £3) £1 + - ^ (£2 - - E 3 ) E 2 +£ kq -^ 

f(E 3] ) + f(E H ) E 3 +£ kq + ^ X 



{E 3 -E l ){E 3 -E 1 ) E 3 +£ ka -^ 



2A 2 E I 2 ( 



1 



kq 



(1) 



k fcq /eg 

The coefficients in the pseudogap dispersion are 

£(k.q)[/(6c-q T )+/(6c-ql)-2] 



kq 

A 

f fry [/oO + + ( k -f )[/<?) - Ml 



!}■ 



(B5) 



- -E 



/(£7 1T ) + /(£7 U ) £ 1+ £ fc9 + ^ 



(-Ei + £ feg - ^)(^2 + £ fe9 - ^)(E 3 + £ kq - ^) (Ei - ^)(Si - £3) ^ + £ fe9 - 



/(E 2T ) + /(E 2i ) £ 2 +£ feg + ^ /(£ 3 t) + /(£3i) + 



kq 



+ 



(£ 2 - E^iE* - E 3 ) E 2 + £ kq -^ (E 3 - Ei ) (E 3 - E\) E 3 + £ kq -^ 

[1 - /(a- qT ) - /(a- q i)] - ^§[1 - /(*© - /(^J)] 



2A 2 



(B6) 



and 



/(6c-q T )+/(£k-ql) 



2M* 



l 2mZ M + £ kq ±*)(E 2 + £ kq ±Z)(E 3 + £ kq - i^) 



+ ^(k-q)(e k -q + M)(- 



1 



1 



— (k-q) + l(6c- q + M) 
m m 



1 



" ^E 1+ £ kq -^ E 2 +£ kq -^ E 3 +£ kq -^ 

£(k • q)(£k-q + M)[/'(£k- qT ) + /'(&-„!)] 



+ 



+ 



+ 



(Ei + 5 fc , - ^)(E 2 + £ kq ^)(E 3 + £ kq 

f(E^ + f(E n ) l ±(^_ tl + ^(E 1 +£ kq + ^) 
(Ei — E 2 )(Ei — E 3 ) 

f(E 2 {) + f(E 2l ) 



(E 2 — Ei)(E 2 — E 3 ) 

f(E 2 i) + f(E 2l ) 
(E 2 — Ei)(E 2 — E 3 ) 

1 



(Ei+£ kq -^) 2 
^ ($k _ q + M )( g2+ g fcg + ka) 

{E 2 +£ kq -^Y 

r (a-q + M)(^3+^ + ^) 



- 6 



6 



Ei+£ kq + ^ 



kq 



E 1 +£ kq - ( 
E 2 +£ kq + ^ 



E 2 + £ kq -^ 
E 3 +£ kq + ^ 



-12^E{[/(&-iT) + /(&-^)] 



(£ 3 + ^-^) 2 

< } (k-q) 2 



3+- 



2) 



3^ , 2[(k- 



E, 



(i) 



+ 



mA 2 

? (i) 2 

^fcq ^kq 



9 J (% + ?ka ) 



+ 



+ -[/(&-«) + /'Kk-„)](k-q) ! 



'( 



3 1- 



kq 

AD 

?kq 



(1) 
fcq 



2£ 



(i) 



inA? 



toAj 

^ 2 ( k-q) 2 }. 



-)) 

kq 



(B7) 
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Here we define u[f = | (l + J^) and v[f = A (l - ||), where - £ kg - k • q/2m. 

* kq ' ^ kn ' 
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